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We study long polymer chains in a poor solvent, confined to the space between two parallel 
hard walls. The walls are energetically neutral and pose only a geometric constraint which changes 
the properties of the coil-globule (or "9-") transition. We find that the 9 temperature increases 
monotonicaUy with the width D between the walls, in contrast to recent claims in the literature. 
Put in a wider context, the problem can be seen as a dimensional cross over in a tricritical point 
of a (f)'^ model. We roughly verify the main scaling properties expected for such a phenomenon, 
but we find also somewhat unexpected very long transients before the asymptotic scaling regions 
are reached. In particular, instead of the expected scaling R ~ N'^^'^ exactly at the (D-dependent) 
theta point we found that R increases less fast than N^^^ , even for extremely long chains. 



I. INTRODUCTION 

Thin films and quasi low dimensional systems are of 
obvious technological interest, in areas ranging from elec- 
tronics to anti-corrosion coatings. Due to this, also crit- 
ical phenomena in systems which are finite in one di- 
rection but infinite in all other ones have been studied 
in quite some detail 0, ■ As the thickness of the film 
decreases, the properties cross over from bulk (3d) be- 
haviour to surface (2d). Actually, in the true thermody- 
namic limit the system should show the scaling typical 
for 2d critical phenomena for all finite thicknesses. But 
since the behaviour must, for finite systems, resemble 
those of 3d critical systems when the thickness is larger 
than the lateral size, the amplitudes in the thermody- 
namic must show some special scaling. All this is known 
as dimensional cross over. The forces exerted by the crit- 
ical fluctuations on the walls are known as Casimir forces 
. The phenomenon exists also for tricritical points, al- 
though it has been studied much less for them. 

Long flexible polymers in very diluted solutions can 
be described as the limit of the 0^ 0{n) vector model 
in the limit n ^ The coil-globule ("theta") tran- 
sition happening as the solvent quality becomes worse 
(typically as temperature is lowered) is in this framework 
described as a tricritical point. Therefore, the problem 
of a polymer confined within the gap between two par- 
allel plane walls is formally described as a dimensional 
cross over, either at a normal critical point (athermal 
polymers) or at a tricritical point (theta polymers), and 
the forces exerted by such a polymer are analogous to 
(critical/tricritical) Casimir forces. 

But polymers have some special features which find 
no close analogy in 0{n) models with n > 0. One of 
them is the fact that the volume occupied by an ather- 
mal polymer of fixed chain length N and confined within 
two athermal walls with distance D does not depend 
monotonously on D If D is larger than the Flory 
radius Rp oc N'^ (with v « 0.5876), then the main effect 
of the confining walls is to reduce the size perpendic- 
ular to the walls, and thus both the gyration and the 
end-to-end radius shrink with decreasing D. But when 



the polymer is strongly compressed {D <C Rp), then the 
main effect is the lateral swelling due to the increased 
excluded volume interaction. In this case, R^ = R\ + i?| 

and the occupied volume R±^ x i?| both increase when D 
is decreased further. 

This observation was the basis for a recent claim 
that the theta temperature Tg of a polymer in a 
poor solvent between two athermal walls should also be 
non-monotonic in D. The argument is essentially that 
the monomer density controls the number of monomer- 
monomer contacts, and thus also Tg: The higher is 
the density, the more effective will be the attractive 
monomer-monomer interaction, and the higher will be 
Tq. Therefore, as D is increased, Tg should start at its 
2d value, increase, go through a maximum, and finally 
decrease in order to reach its 3d value when Z) — > oo. 
This predicted behaviour was then supported by exact 
enumeration studies of short chains {N < 20) on the 
simple cubic lattice (see Fig. 6 below). 

In the present paper we show by means of rather ex- 
tensive Monte Carlo simulations that this is not correct, 
and that Tg increases monotonicaUy with D. This is es- 
sentially what one would have expected from a dimen- 
sional cross over of a tricritical point. On the other 
hand, some of the details of this cross over are somewhat 
surprising. In particular, we find extremely long tran- 
sients. At T = Tg{D), the asymptotic scaling R\\ ^ N""-^ 
with j/e_2 = 4/7 being the Flory exponent for 2d theta 
polymers 0, IBl is not seen even when Rp/D w 100, 
although we have no reason to doubt that it will hold 
for Rp/D oo. The same is true for the free energy: 
Although we have no doubt that the scaling appropri- 
ate for 2d theta polymers will apply asymptotically at 
T — Tg{D), it is not yet seen in the simulations. 

These simulations are done with the pruned-enriched 
Rosenbluth method (PERM) 9] which is ideally suited 
for this purpose. It allows to study extremely long chains 
(for D ~ 60 we went up to iV = 600, 000) with very high 
statistics, and it gives immediately very precise estimates 
of free energies. Throughout the paper we shall model 
the polymers by self avoiding walks on the simple cubic 
lattice with attractive energy — e between non-bonded 
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neighbouring monomers. Instead of quoting tempera- 
tures or values of e, we shall describe thermal effects in 
terms of the Boltzmann factor 



(1) 



per contact. The partition sum is therefore 

ZAr((J,i?) = ^CAr^„,(i?)q™ (2) 

m 

where CN,m{D) is the number of walks with N steps and 
m monomer-monomer contacts. The width D is defined 
such that D = 1 corresponds to the standard square lat- 
tice. For large D we used hashing as described e.g. in 
|10| in order to minimize storage demands. 

II. NUMERICAL RESULTS 

In the following, we will denote by qe{D) the value of 
the Boltzmann factor at the true quasi-2d theta point. 
The Boltzmann factor for 3-d theta polymers in the bulk 
is then = lim/j^+oo qe{D), while the Boltzmann factor 
for strictly 2-d theta polymers on the square lattice is 
Qg^^ = qs{D = 1). The same notation is used for theta 
temperatures and for growth constants (inverse critical 
fugacities). Growth constants at temperatures different 
from the collapse point will be denoted as fi{q,D), so 

that e.g. = lim£)^oo K^f\D). 

Exactly at the 3-d tricritical (theta) point, we can as- 
sume finite size scaling ansatzes for the partition sum 
and for the rms. end-to-end distances both parallel and 
perpendicular to the walls: 



ZN{qf\D)^[^ifYHD|RF{N)), 



(3) 



{{xN - xnf + (yjv - yof) 
RFiN)^\liD/RpiN)) (4) 



and 

(5) 

up to logarithmic corrections P. ITolllllll^ . Here, Rf{N) 
is the Flory radius (rms. end-to-end distance of chains 
in the bulk) which scales at the theta point like N^^^, 
again up to logarithmic corrections. The scaling func- 
tions ^(z), (z), and '^±{z) are finite and non-zero in 

(3) 

the limit z — > oo. Finally, the growth constant /i^ is 
a non-universal constant which for the present model is 
5.0479050 ± 0.0000005 -h(g^^'- 1.3087) x 1.616 [HHIl. 

If the theta temperature T0{D) increases with 1/D, 
sufficiently long polymers a,t q — qg ' will be collapsed 
for finite D (i.e., i?A,(gf £>) - N'^/^), while they will 
be swollen (i?jv('7f \ £>) ^ N^^^) iiTe{D) decreases with 
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FIG. 1: Squared end-to-end distances parallel to the walls, 
divided by N^^'^ , plotted against A''. For self avoiding walks 
in 2d one has ~ N^^^, i.e. the walks are swollen coils if 
the lines become horizontal for large A'^. 



1/D. In Fig. 1 we plot R%^^^{q'^g^\ D)/N^/'^ against N, 
for various values of D. We see that all curves become 
horizontal for large N, i.e. all chains are swollen. This 
is in contradiction to the claim of @. The reason why 
the heuristics leading to this claim were wrong is indeed 
quite clear: The location of the theta point is determined 
by polymers whose Flory radius is much larger than D, 
while the anomaly noticed by van Vliet et al. _5j concerns 
only polymers with Rp ~ D. 

From Eq.lQ we see that \l/||(z) ~ z~^/^ for z ^ 0, 

in order to obtain R]^ \\(qf \ D) ~ A(D)N'^^'^ with some 
amplitude A{D). The ZJ-dependence of the amplitude is 
then also fixed by Eq.Q), 



(6) 



To check this, we plot in Fig. 2 the amplitudes obtained 
by fitting horizontal lines to the large- TV data in Fig. 1. 
We see that the behaviour is roughly as predicted in 
Eq.®, but not quite. There are obviously substantial 
logarithmic corrections (similarly large corrections are 
also seen in theta polymers in the bulk, see [ol fl^ [T3|). 
In the eyeball fit shown in Fig. 2 these are described by 
a factor cx ln(_D)^/*, but this is just done to guide the 
eyes. A similarly good fit would have been obtained with 
a pure power law A{D) ~ D^^-^'^. 

For all values of D, the partition sum for N ^ £)i/''2 jg 
also compatible with the scaling behaviour expected for 
2-d SAWs (i^2 = 3/4 is the Flory exponent in 2 dimen- 
sions), 



Zj,{ql,'\D)^ZoiD)fl{ql,'\DfN'<^ 



(7) 



where 72 — 43/32 = 1.34375 and where 
^ksTg In ji{qg ,D) is the free energy per monomer 
in the thermodynamic limit. This is compatible with 
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FIG. 2: Amplitudes for the end-to-end distance parallel to 
the walls (see Eq.Q) plotted against D. The smooth line 
is an eyeball fit with a power law oc D~^^^ modified by a 
logarithmic term. 



FIG. 3: Dependence of the effective connectivity constant 
fj,{q,D) on D, for fixed q = qg^'' . The straight line gives the 
asymptotic behaviour ~ which is obviously modified 

by logarithmic corrections. Error bars on the MC data are 
much smaller than the symbol sizes. 



Eq. ( |3Jl , provided that 

(3) 

Me 



1-^ 



Zo{D) - 1)2-272 



and 



<i>(z) exp(-o/z^) 



for 







(8) 
(9) 

(10) 



Eq. © can also be derived from field theory Q, |^ |^ ^3 j 
using De Gennes' mapping onto the 0{N = 0) model. 
Numerical values for /ig^"* — ^{qg'^\ D) are shown in Fig. 3. 
They are obviously in agreement with Eq. |(SJ| for large 
values of -D, although there are large corrections for finite 
D. 

Up to now we have only discussed the behaviour at the 
theta point of the free {D = oo) polymer. Let us finally 
study the behaviour at and in the vicinity of the theta 
points for the quasi-2d systems at fixed finite D. Accord- 
ing to the above results, this corresponds, for any finite 
D, to temperatures where the free 3-d polymer would 
be collapsed. Asymptotically, for iV ^ oo, the behaviour 
exactly at the tricritical point should be that for 2-d theta 
polymers, 



ZN{qe{D),D)^[ii{qe{D),D)]''N'" 



and 



RN,\\{qe{D),D) 



(11) 



(12) 



Here, i'g^2 — 4/7 and 761^2 = 8/7 are the Flory and en- 
tropic exponents for 2-d theta polymers. 

In Figs. 4 and 5 we show R% ^^ {q,D)/N^/'^ versus N, 
for several values of q close to (tri-)criticality, and for 



D — 5 (Fig. 4) resp. D — 60 (Fig. 5). According to 
Eq. (|12() , one of these curves in each graph should become 
horizontal for iV — > 00. Naively, one might expect this 
regime to set in when the (3-d) Flory radius is roughly 
equal to D, i.e. for iV w 25 in Fig. 4 and for N w 4000 
in Fig. 5. We do not see this, although our values of 
N are much larger than these, by factors more than one 
hundred. Although we cannot pin down precisely qe{D) 
due to this, we definitely see that Rn increases at the (D- 
dependent) theta point - and for numerically accessible 
values of - slower than \/N (straight lines in Figs. 4 
and 5). Thus instead of being swollen, as predicted by 
theory, the polymers seem to be collapsed at the 3-d theta 
point. Similar results were found for all other values of 
D, and they are corroborated by the results obtained 
for Zjv- Also for Z^, the asymptotic behaviour stated 
in Ea. ljllll was not seen for any of the chains (data not 
shown) . 

While these results might look very surprising at first 
sight, it is indeed not too difficult to understand them 
heuristically. As we said, free 3-d chains would be slightly 
collapsed at the temperatures shown in Figs. 4 and 5. In 
a blob picture, these chains are therefore chains of blobs, 
each one of size « D, and each representing a short 3-d 
polymer slightly below the theta point. It is well known 
that the 3-d theta collapse of finite chains hap- 
pens at an effective ("Boyle") temperature Tg{N) which 
is lower than the true Tg , and that concatenating such 
chains at Te{N) leads to a longer chain which is much 
more collapsed than its short constituents: For tempcr- 

(3) 

atures slightly below Tg , there is a regime in N where 
dRpf/dN < 0, i.e. chains actually shrink when more 
monomers are attached to them! For the present problem 
this means that each blob might be swollen, but when two 
blobs are brought in contact, they might not repel each 
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FIG. 4: Average squared end-to-end distances parallel to the 
walls for nearly tricritical polymers in a slit of width D = 
5. More precisely, R%/N^^'' is plotted against InA'^, which 
would lead to a horizontal curve if Ea. ll2ll were satisfied for 
all A^. The straight line indicates, in contrast, a non-swollen 
behaviour R% oc A'^. 



FIG. 6: Lower limit estimates of qg{D), plotted against 1/D. 
The straight line is q = 1.305 -I- 0.53/Z). The insert shows 
same of these data plotted against D, together with the data 
from Ref. [6] for comparison. 
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FIG. 5: Same as Fig. H but for D = 60. Notice that the 
longest chains in this case have N = 600, 000. 



other (as in a swollen polymer), but they rather attract 
each other. Obviously this is what happens. For much 
longer chains than those we can simulate, this would fi- 
nally stop: When too many blobs are penetrating each 
other, repulsion finally dominates again and the chains 
behave as ordinary 2-d theta polymers. 

According to this discussion we can give only lower 
bounds on qe{D). They are given by those values of 
q for which the curves in plots like Figs. 4 and 5 be- 
come horizontal. These bounds, which should however 
be not too far from the true values of qg(D), are shown 
in Fig. 6. We see an essentially linear increase with 
1/D. Extrapolating linearly to l/D ^ we find a value 
q = 1.305 ±0.001 which is close to, but definitely smaller 
than, ^ = 1.3087 ± 0.0003 0. This confirms that the 



curves in Figs. 4 and 5 which become horizontal at the 
largest values of A^ are not yet the critical ones, and that 
the true asymptotic behaviour is not yet seen in these 
figures. 



III. DISCUSSION 

On the one hand we have shown that some recent 
claims about re-entrant behaviour of polymer collapse in 
restricted geometries are wrong. The theta collapse of a 
polymer confined to the space between two parallel walls 
is very much as expected from a tricritical behaviour in a 
film (quasi-2d) geometry. In particular, the collapse tem- 
perature is, for any distance D between the walls, shifted 
to temperatures lower than the theta temperature of free 
3-d polymers, and the decrease of the free energy with 
decreasing D is as predicted from the theory of tricritical 
Casimir effects. 

But on the other hand, we found that the detailed be- 
haviour at the true quasi-2d theta collapse is - for any 
chain length we could simulate, and presumably also for 
any chain length realistic in any foreseeable experiment 
- rather different from the predicted one. This is related 
to the fact that the upper critical dimension of theta col- 
lapse (as of any other tricritical phenomenon) is d = 3. 
Therefore long chains are, at the collapse point, com- 
posed of blobs which are essentially free random walks 
but which slightly attract each other. 

In view of the analogy between theta collapse and 
other tricritical phenomena, it is of interest to specu- 
late whether similar anomalies should be expected also 
for the latter. For tricritical Ising or Potts models one 
might then expect that the effective correlation length 
exponent is not given by the true tricritical exponent for 
d = 2 (which is larger than 1/2 for all these models), but 
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has a value < 1/2. To our knowledge this has neither 
been predicted so far, nor has it been seen in simulations 
or in real experiments. 
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